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We propose a digital quantum simulator of non-Abelian pure-gauge models with a superconduct¬ 
ing circuit setup. Within the framework of quantum link models, we build a minimal instance of 
a pure SU (2) gauge theory, using triangular plaquettes involving geometric frustration. This real¬ 
ization is the least demanding, in terms of quantum simulation resources, of a non-Abelian gauge 
dynamics. We present two superconducting architectures that can host the quantum simulation, 
estimating the requirements needed to run possible experiments. The proposal establishes a path 
to the experimental simulation of non-Abelian physics with solid-state quantum platforms. 

PACS numbers: 03.67.Ac, 03.67.Lx, 85.25.Cp,11.15.Ha 


Gauge invariance is a central concept in modern 
physics, being at the core of the standard model of ele¬ 
mentary particle physics. In particular, invariances with 
respect to SU{2) and SU{3) gauge symmetries charac¬ 
terize the weak interaction and quantum chromodynam¬ 
ics ID HI- In this sense, gauge theories represent a cor¬ 
nerstone in our understanding of the physical world and 
lie at the heart of diverse phenomena, such as the quark- 
gluon plasma or quantum spin liquids. In condensed mat¬ 
ter physics, SU (2) gauge fields can also emerge dynam¬ 
ically in relation to exotic many-body phenomena, like 
quantum Hall systems, frustrated magnets, or supercon¬ 
ductors [sHg. 

Lattice gauge theories (LGT) are non-perturbative dis¬ 
crete formulations that contribute to the analysis of key 
features of these models, such as color confinement or 
chiral symmetry breaking. Starting from the seminal 
work by Wilson in 1974 EUD, LGT have attracted a 
significant attention across several branches of theoreti¬ 
cal physics. In the last decades [8], quantum Montecarlo 
simulations have achieved unprecedented accuracies in 
determining the whole hadronic spectrum of the standard 
model. However, understanding its full phase diagram 
from first principles, or simulating dynamical processes, 
remain out of reach of current numerical computations. 

Quantum simulators [9] provide a new approach 
to solve complex long-standing problems in quantum 
physics. In a quantum simulator, a proper encoding 
of LGT can allow for the retrieval of information about 
ground state and dynamics, in a wide range of regimes. 
Previous works have considered Abelian [T0]-[22] and non- 
Abelian LGT in optical lattices [23ll25] (see also EHEi] 
and references therein), both as analog and digital sim¬ 
ulations [28]. In these implementations, matter-gauge 
interactions are modeled as a second-order process from 
Hubbard-like interactions, while the simulation of pure- 
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FIG. 1. (Color online) a) Six tunable-coupling transmon 
qubits coupled to a single microwave resonator, b) Six Xmon 
qubits on a triangular geometry, coupled to a central one. The 
box 1 in the scheme is implicitly repeated for the sides 2 and 3. 
Both setups can encode the dynamics of the SU (2) triangular 
plaquette model schematized in c), where the left and right 
gauge degrees of freedom are explicitly depicted. 


gauge dynamics remains more demanding. 

In the last years, superconducting circuits have proven 
to be reliable devices that can host quantum information 
and simulation processes m- The possibility to perform 
quantum gates with high fidelities, together with high 
coherence times, makes them ideal devices for the real¬ 
ization of digital quantum simulations [3Qll34] , previously 
considered in ion-trap systems |35H37] . 

In this Letter, we propose a digital quantum simula¬ 
tion of a non-Abelian dynamical SU (2) gauge theory in 
a superconducting device. We start by building a mini¬ 
mal setup, based on a triangular lattice, that can encode 
pure-gauge dynamics. The degrees of freedom of a sin¬ 
gle triangular plaquette of this lattice are encoded into 
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qubits. We propose two implementations of this quantum 
simulator, using two different superconducting circuit ar¬ 
chitectures, as depicted in Fig.[2 We consider a setup in 
which six tunable-coupling transmon qubits are coupled 
to a single microwave resonator, and a device where six 
capacitively-coupled Xmon qubits stand on a triangular 
geometry, coupled to a central auxiliary one. We com¬ 
pute experimental requirements necessary to perform the 
simulation on one plaquette, and provide arguments for 
scaling to large lattices. 

Lattice gauge theories.- LGT are discretized versions 
of a gauge theory. In a conventional approach of lat¬ 
tice gauge theories, space-time is discretized while ensur¬ 
ing a covariant formulation of the theory. In quantum 
simulations, there is no direct access to the time direc¬ 
tion, that is fixed and continuous. Hence, the equivalent 
Hamiltonian formulation of lattice models is used. In 
this case, the space is discretized and the action of a 
local gauge invariant Hamiltonian characterizes a contin¬ 
uous dynamical evolution. This fictitious asymmetry of 
the time direction forces us to define a physical Hilbert 
space or Gauss law of the model. We make use of a 
set of discretized space points x, and SU{N) operators 
u{x^ jl) associated to each link, connecting two adjacent 
sites (x, x-h/i). A single u{x^ jl) operator has two implicit 
color indices ((a,/d), Ucxjs^ that connect the a-th compo¬ 
nent of the fermionic field at position x with the P-th 
component of the field at position x + /i. 

In order to construct a non-Abelian SU{N) theory 
built out of these u{x^jl) operators, one has to build 
a Hamiltonian that is invariant under a generic gauge 
transformation 

u{x,r) ( 1 ) 

where r“ are the — 1) generators of the SU{N) al- 
gebra in the corresponding representation, and 0^ the 
relative phase angles. It is straightforward to verify that 
the following Hamiltonian on a square lattice is gauge 
invariant 

H = J'^^Tt[u{x^ jl)u{x-\-jl^i>)u\x-\-i>, jl)u\x^ z>)]+H.c., 

X 

( 2 ) 

where the trace is performed over the color indices, J has 
energy units, and the directions jl and z> span the two- 
dimensional lattice. This Hamiltonian is a pure gauge 
operator, as it does not involve any fermionic operator 
and it is the minimum instance of a Wilson loop around 

a plaquette mm- 

We focus now on the construction of a non-Abelian 
SU{2) quantum link model that mimics this gauge- 
invariant behavior in a finite-dimensional Hilbert space, 
suitable for quantum simulations. Notice that the frame¬ 
work of the quantum link models is valid for any com¬ 
pact Lie group [38]. We use a four-dimensional Hilbert 
space or, equivalently, two qubits associated with each 
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FIG. 2. (Color online) a), b), and c) Different gauge invariant 
sectors for a triangular plaquette, together with the sector 
degeneracy and the action of the Hamiltonian Ht upon them, 
d) Extended lattice with periodic boundary conditions. This 
lattice allows for the G‘^{x) = 0 sector. 


link. One can define link operators U{x^jl) acting on 
this four-dimensional Hilbert space with the associated 
gauge generators G^{x) that satisfy the SU{2) algebra 

[G“ {x ), (y)] = i8s;g E] eabcG'" {x). (3) 


To define the set of operators, we use a quantum link 
formulation [38] , where the Hilbert space of a link is given 
by a set of bosonic modes Caj {x, jl) that implement the 
Schwinger representation of the SU (2) algebra, with a G 
{t, i}, acting on two different sites j G {L^R} on the 
link between the two adjacent vertices x and x 1- ji. In 
this space, we build two sets of right and left generators 
I/^, on this finite-dimensional Hilbert space, using 
the bosonic modes c^j (^, jl) [39] . We define the gauge 
generators as G^(;r) = L^{x, t>){x — i), z>), where 

the sum over |z^| is taken among all the links of the lattice 
converging to a single vertex x. 

The representations of SU(2) are quasi-real, therefore, 
ordinary and their dual representation (i.e., particle and 
antiparticle) are equivalent. Therefore, there are two 
multiplets with well-defined SU{2) transformations [39] . 
One can finally define the link operators, acting on the 
finite-dimensional Hilbert space of one link, in terms of 
Schwinger bosons on a given link [47]-[49], 


Ua.g — 



ie\L J 


-ic^Rj 


(4) 


These operators U have the following commutation rules 
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with the left and right operators 

(7^ 

[U{y,v),R^{x,ij)] = -U{y,fi)—6s;ySofi, 

' ( 5 ) 

[U{y,v),L‘^{x,fi)] = —U{y,jl)5sg5{,(„ 

where we have defined the usual Pauli matrices cr^, see 
also [39], while is defined as the identity operator. The 
commutation rules for a general gauge transformation 
follow in a straightforward way, 



JJ e-iS‘'iy)G‘^iv)U{x, fl) JJ 

y Z 

= (6) 

where the products over y and z are extended over the 
whole lattice. This equation shows that quantum link 
models are formulations of gauge-invariant models. In 
fact, one can notice how the action of a gauge transfor¬ 
mation in Eq. mimics Eq. 0- 

A generic state in the local Hilbert space 
^ quantum link is given by 
the occupation of the operators R{l)' Since the 
total occupation per link, + R^r^ is 

a constant of motion [see Eq. Q], we restrict to the 
subspace with one occupied mode per link. 

Notice that matter-gauge interactions in 1 + 1 dimen¬ 
sions can be derived [39] in a second-order perturbation 
theory, considering extensions of previous works [50] [5l]. 
In the following, we will focus instead on pure-gauge two- 
dimensional interactions. 

Triangular lattice.- The continuum limit of quantum 
link models and their thermodynamical properties are 
topics under active research. In related condensed matter 
models, the importance of the lattice geometry has been 
shown in fundamental aspect of the phase diagram, as 
the existence of confinement and de-confinement phases 
[5^ l53] . Quantum simulations on large lattices may pro¬ 
vide insights on the continuum limit of LGT, and its 
dependence on the geometry of the underlying lattice. 
In this spirit, we consider a minimal implementation of a 
pure SU (2) invariant model in a triangular lattice, by us¬ 
ing triangular plaquettes, as depicted in Eig. [^. In this 
case, the pure-gauge Hamiltonian on a single plaquette 
reads 


Ht = —J Tr \U (x, jT)U{x + /i, 0)17 (x + /i -1- z>, —ja — z>)]. 

(7) 

This interaction corresponds to the magnetic term of a 
gauge invariant dynamics, which acts on closed loops. 
We focus on the magnetic term, since the representation 
of the electric counterpart is trivial. 

Since this Hamiltonian commutes with gauge genera¬ 
tors, = 0, Vx, a, these are constant of mo¬ 

tion, and one can define different gauge sectors with the 
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FIG. 3. (Color online) The relative gauge deviation of the 
digitally simulated gauge value {{)d), versus the ideal value 

((>/), E = {{^,GHx)}r - {E,GHx)}o) as 

a function of the digital steps N and the simulated phase 
4> = Jt. The error decreases with large A, and small 0. 
The contour plot is interpolated from numerical data at A = 
{1, 2, 3,4}. The initial state is chosen as depicted in Fig. [^. 


initial values of the generators. Due to the triangular ge¬ 
ometry and having one occupied mode per link, there are 
no global states that have zero eigenvalue for the gauge 
generators G^{x) at every vertex in a single plaquette. 
In general, the absence of the zero-eigenvalue sector de¬ 
pends on the topology of the lattice, and is avoided, for 
example, in the periodic lattice of Fig. Ht- The different 
gauge sectors can be characterized by the eigenvalues of 
the operator in Fig.[^, 

[^, and 0, see the supplemental material for additional 
details [39] . 

The local Hilbert space of a link is four dimensional, 
and it can be faithfully spanned by two qubits, called 
“position” (jpQg and “spin” qubit cr^. In this subspace, 
it is useful to define the operators F^ = ^pos^m? = 
^pos^m 5 such that the total Hamiltonian is written as [39] 


Ht = 


-^{r? 


pO pO 
12^- 23J- 31 




(8) 


abc 


-E 


pO pa pa 
J- 12^- 23J- 31 


r a pO pa I pa pa pO 
12-'- 23-*- 31 I 2.“^-'- 


12-L 23-l 31 


Superconducting circuit implementation.- In order to 
simulate the interaction of Eq. one can decompose 
its dynamics in terms of many-body monomials, and im¬ 
plement them sequentially with a digitized approxima¬ 
tion [54]. In a digital approach, one decomposes the dy¬ 
namics of a Hamiltonian H = by implement¬ 
ing its components stepwise, ~ 

(here and in the following h = 1), for a total of 
m X N gates, with an approximation error that goes to 
zero as the number of repetitions A grows. In a practical 
experiment, each quantum gate will be affected by 

a given error e/^. By piling up sequences of such gates, 
for small gate errors <C 1, the total protocol will be af- 
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fected by a global error, which is approximately the sum 
e ^ e/c- This model for error accumulation has been 
proved in recent experiments [32l |33] . 

The effects of digitization in the simulation of LGT 
can be observed in Fig. We have numerically inte¬ 
grated a Schrodinger equation regulated by the Hamil¬ 
tonian in Eq. (§, choosing as initial state one of the 
12 states shown in Fig. [^, for different simulated 
phases (p = Jt [39]. We compute both the evolu¬ 
tion with the ideal Hamiltonian and with a digital se¬ 
quence, in which we act with a single monomial at 
a time, repeating the protocol N times. As a fig¬ 
ure of merit, we compute the relative deviation E = 

{{E,GHx))i - {E,GHx))d) /{E,GHx))i from the 
ideal value of the gauge invariant, where {)i{d) stand 
for average values computed on the ideal (digitally simu¬ 
lated) state. The deviation goes to zero for large N and 
small phases, defining surfaces at given error tolerance in 
the N-cf) space. 

To simulate the pure-gauge interaction in a single tri¬ 
angular plaquette, we first consider a setup in which six 
tunable-coupling transmon qubits are coupled to a sin¬ 
gle microwave resonator [55l [56] . Each tunable-coupling 
qubit is built using three superconducting islands, con¬ 
nected by two SQUID loops. Acting on these loops 
with magnetic fluxes, one can modify the coupling of 
the qubits with the resonator, without changing their 
transition frequencies. For further clarifications of the 
experimental setup involving tunable-coupling transmon 
qubits, see the Supplemental Material [39]. By threading 
with magnetic fluxes at high frequencies, one can drive 
simultaneous red and blue detuned sidebands, and per¬ 
form collective gates m- Each many-body operator can 
be realized as a sequence of collective and single-qubit 
gates, 

Un = E.<,- g-i(^/4) ^ (9) 

where the indices i and j run from 1 to A". In this way, 
one can implement a generic (up to 2N local rotations) 
A-body evolution operator Un = - with 

{m,n, .../c} G {x,^,z}, using 2 collective gates, and a 
number of single qubit gates which is upper bounded by 
2A + 1, counting the single qubit gate in Eq. and the 
2A rotations necessary to map Uw to any A-body oper¬ 
ator. The simulation for one digital step will amount to 
implementing 32 collective gates and a number of single 
qubit gates which is upper bounded by 184 [39] . 

We consider now the architecture of Fig. [^, in which 
six Xmon qubits in a triangular geometry are capacitively 
coupled with an additional central ancillary qubit [58] . 
In this case, the collective interactions Un can be de¬ 
composed and performed with pairwise C-phase gates, 
using the central ancillary qubit to mediate non-nearest 
interactions. In this way, the quantum simulation of one 
digital step of the Hamiltonian in Eq. ^ will amount to 



FIG. 4. (Color online) Digital quantum simulation of pure 
gauge dynamics for a) 2 and b) 3 digital steps. The initial 
state I To) is chosen as in Fig. [^. Periodic oscillations be¬ 
tween the position and qubit degrees of freedom at the ver¬ 
tices are witnessed by |(T/|To)|^. The theoretical fidelity 
|(T/|Ti:))|^ increases with the number of digital steps, to¬ 
gether with the width of the fidelity bands due to experi¬ 
mental imperfections. 


realize 168 C-Phase gates and a number of single-qubit 
rotations which is upper bounded by 520. For addi¬ 
tional details on the experimental setup and gate counts, 
see [39] . 

To estimate the effects of experimental imprefections 
on the simulations, we plot in Fig. [^the dynamics of a 
single plaquette, starting from a = 9/4 state. The 

observed periodic dynamical oscillations between the two 
paired states, depicted in Fig. It, are signaled by the 
periodic behavior of the overlap |(T/|To)p. The simu¬ 
lation is run for A = 2 and A = 3 steps, necessary to 
observe one flip between the two states and one full oscil¬ 
lation, respectively. On top of these oscillations, we plot 
fidelity bands that estimate the fidelity cap due to accu¬ 
mulated gate errors in the digital protocol, in the case 
of the two presented setups in Fig.[^ with collective and 
two-qubit gates. In order to see coherent oscillations, the 
fidelity of the two-body gates has to be one order of mag¬ 
nitude better with respect to the collective ones, with an 
error window of approximately 0.05 — 0.01% in the col¬ 
lective gate case and a 0.005 — 0.001% for the G-phase 
gate setup. We assume a factor 1/20 for the error loss 
on a single qubit gate. From the plots, it can be seen 
that, while the digital fidelity |(T/|Ti:))p increases, the 
fidelity bands due to experimental errors broaden, affect¬ 
ing the simulation. The intersection of the digital fidelity 
line with the bands marks a regime dominated by exper¬ 
imental imperfections, for small 0, and one in which the 
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error of the digital expansion prevails, for large (j). 

The proposed simulation can be extended to larger lat¬ 
tices, where one could analyze dynamical properties of 
the confinement-deconfinement phase transition, which 
can be traced back to the breaking of the global center 
symmetry m in groups with a non-trivial center group, 
like the SU (2) case. Moreover, such a LGT simulator can 
be used to perform quench experiments in a ladder con¬ 
figuration, analyzing breaking of gluon strings between 
pairs of particles and anti-particles [6Qh6^ . 

In order to scale the quantum simulation to large qubit 
lattices, one has to consider that the accumulated gate 
error does not depend on the size of the lattice [32j |33] , 
and only on the number of gates. The quantum resources 
necessary to run the simulation scale polynomially in the 
lattice size, and sub-polynomially in the digital error [63] . 
making the whole protocol efficient. With a current to¬ 
tal number of gates that exceeds 1000 [32l|3lj, we believe 
that in the near future simulations of LGT on large scales 
will be feasible. Furthermore, the possibility of perform¬ 
ing quantum error correction for digital quantum simu¬ 
lations may drastically increase its effectiveness [64] . 
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SUPPLEMENTAL MATERIAL 

In this supplemental material, we perform additional analysis that support the results obtained in the main article. 


Explicit SU{2) gauge operators in terms of Schwinger modes 

In this section, we explicitly derive the link operators in terms of Schwinger bosons. Following [T], we define the 
operators as right and left generators, acting on the finite-dimensional Hilbert space of the link. 


0 "^ 

{a,/3}G{t A} 

cYx,fl)-^Ci3L{x,ll), 

{a,/3}G{t ,4} 

where c^j are bosons that implement the Schwinger representation of the SU (2) algebra, and the Pauli matrices 


( 10 ) 




1 0 
0 1 


; = 


0 1 
1 0 


; = 


0 -i 
i 0 


; = 


1 0 

0 -1 


( 11 ) 


The bosonic operators act on two different sites p G {L^R} on the link between the two adjacent sites x and x P fi 
and obey the SU (2) commutation rules 

R^]=iY [-R“, L^] = 0, [L“, L^] =i Y 


( 12 ) 
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taking into account that they commute on different links (i.e., different x, fi). With the Schwinger representation, 
there are two multiplets with well-defined SU (2) commutation relations 






y^.'^ad^BL 




E 




B 


J DO 
"ai?’ ^ 


Z^S-R 9 ’ 




■E 

7 V /3 


Given these two multiplets, we can define the link operator as 

UaB = CaLC^ji + ^ 


(13) 


(14) 


We stress that these quantum link operators, or gauge fields, are associated with links between lattice points, 
transporting color information between the two points they connect. The two indices of the link operator Uap{x,fi) 
are identified with the two ends of the link (x, fi). The left (right) index is associated with the beginning (end) of the 
link as shown in Fig. The gauge transformation acts on U{x^ jl) according to 

U{x,ii) = (15) 


Since in general 0^{x) and 0^{x + ft) are different, our gauge-invariant equations will require invariance under right 
multiplication and left multiplication separately. An arbitrary gauge transformation can be built from individual 
gauge transformations at the point of the lattice. Therefore, it is sufficient to consider just a gauge transformation at 
the lattice site x. Every link emanating from the point is affected by the gauge transformation generated by 

(x) = ^1^^ {x - jl,jl) + {x,jl), (16) 

A 


and the total color (x, jl) carried by a link is just the difference: (x, jl) — (x, jl). 

The sum over non-Abelian electric fields emanating from a single point of the lattice is the lattice analog of the 
electromagnetic divergence of the electric field, V • E^. However, because the electric field itself varies from (x^ji) 
to (x, jl) along the link, there is an additional contribution given by the total color carried by the link. Thus, the 
generator is recast into: (x) = V • (x) — ^ (x, jl). The gauge invariance of the physical sector is defined 

by (x) |T) = 0, i.e., at every vertex the state transforms locally as a singlet of SU{2). If we identify the local 
color density {x^ jl) = ^ {x^ jl), then V • E^ (x) = {x, jl) is the usual Gauss’ law. These are properties of 

non-Abelian lattice gauge theories, independent of their representation via quantum link models. 


Local Hilbert space of SU{2) quantum links 

The local degrees of freedom in a SU{2) link are encoded in the quantum state of four different bosonic modes 
^ . Different representations of this local Hilbert space are given by the occupation of these modes. 

The Hilbert space with one excitation in total is four-dimensional, and it can be decomposed in a spin degree of 
freedom with the SU (2) symmetry group and a position degree of freedom. 

A SU (2) gauge invariant interaction commutes with the local generator of the SU (2) at any vertex. As a conse¬ 
quence, the total angular momentum at any vertex of the lattice is a conserved quantity of the dynamics. In the 
simplest example of two quantum links, i.e. just one vertex, there are only two possible gauge sectors, the singlet 
sector where the total angular momentum around the vertex is zero, and the triplet sector where the total angular 
momentum is equal to one. 


Reduction to qubits - Triangular lattice 

In the main text, we have considered the local four-dimensional Hilbert space spanned by the occupation of four 
different modes \n^L^ '^ir)- By choosing the subspace with one excitation per link, we can faithfully represent 








Local Hilbert space of a SU(2) quantum link 
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FIG. 5. (Color online) The four dimensional Hilbert space of the minimal representation of a SU{2) quantum link. Below, the 
two different gauge sectors in the Hilbert space of two quantum links, i.e. the singlet and the triplet sector. 


its quantum state with two qubits, the position and spin qubits, |pos, m^). The excitations of these qubits are related 
to the left and right modes through pos = {n^R + ti^r) — {n^R + + %l) — + ^t^)- 

With these definitions, the left and right operators defined in the main text are recast into 


= 


^al3 


/ . ^aR~ 


C^R = 


7^ + (7^ (7^ 

^pos ^ ^pos ^rn 


TO- - W — 

^ ~ Z_^ ^aL 9 — 


(7 — (7 rr 

"^pos "^pos 


{a,/3}G{t ,t} 


2 ’ 


JJ = I ] 

y^lLC^R — C^l^IR ^ILC^r + C^l^^R J 


I^C^rcIr + C^L^iR cL 


: 0_0 


rv 




^pos^m ^^pos^m 


(jX _o — j(jy (j^ 

^pos^m ‘'^pos^m 


(17) 


with = cTpQgCrJI^, = cr^os^^m- The Hamiltonian for one plaquette can then be written down as 


Hr = — JTr [U{x, jl)U{x + /i, t))U{x + /i + z>, —jl — z>)] . 


(18) 


In order to identify all the terms of the trace one can use the condition tr[r‘^r^T'^] ^ 0. The condition is met if and 
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only if 

a = /3 = 7 = 0, 

a 7^ /3 7^ 7 7^ 0, 

a = 0; = 7 7^ 0, 

= 0; a = 7 7^ 0, 

7 = 0;/3 = a7^0. (19) 

The first condition gives 1 term, the second 6 terms, and 9 terms for the last three equations, for a total of 16 terms 


IV [!7i2C/23!73i] = 1 1 r? 2 r 03 r 0 i + ^ ^ [r?2r^3r3“i + r“ 2 r 03 r 3 “i + \, ( 20 ) 

V abc a ) 

where we have used the subscripts 12, 23, 31 to represent the links (x, /i), (x + /i, z>), (x + /i + z>, —fl — z>), respectively. 
The pure-gauge Hamiltonian can then be written in terms of T matrices as 


^ = -^\ r?273r^i + E^«f'-r“273rii - E [r? 

V abc a 


r a -pa I pa pO pa 1 pa pa p' 
23-*- 31 I 12-'- 23-'- 31 1 12-'- 23-'- , 


12-L 23-L 31 


12-L 23-L 31 


( 21 ) 


Then, one can explicitly write this Hamiltonian in terms of up to six-body interactions between qubits, for instance 

^0 pO pO 


pU pU pU _ X 

12^- 23 J- 31 —<^posl2^pos23^pos31 


pa p6 pc _ y y y 

J- 19J- oa-L “^1 — 


2^p 


l^p 


12 ^ 23^ 31 —^posl2^pos23^pos31^ml2^m23^m31 

E pO pa pa _ X y y \ ^ a a 

J- I2J-23J-31 —^posl2%os23^pos31 / ,^m23^m31- 


( 22 ) 


Experimental setups 


In this section we provide a brief overview of the experimental setups considered for the implementation of the 
proposed simulator. The first setup considered in the main text is comprised of six transmon qubits provided with 
tunable coupling to a single microwave resonator. A transmon qubit is a special instance of a Cooper pair box [2], 
operating at specific regimes. 

A Cooper pair box is made of a piece of a superconducting island connected to a voltage source through a capacitor, 
and grounded via a Josephson junction. The temperature of these superconducting systems is low enough so that 
the quantum state of the setup can be defined by the number of Cooper pairs on the superconducting island. In 
order to add tunability and control to this quantum system, the single Josephson junction of the Cooper pair box 
is substituted with a SQUID loop, threaed with magnetic fluxes. These superconducting systems can be coupled 
to superconducting microwave resonators, which can be used to induce, via electromagnetic interactions, tunnelling 
of Cooper pairs in and out the superconducting island, as well as reading out the state of the system [3l |4]. The 
Hamiltonian of the setup comprising a Cooper pair box coupled to a microwave resonator reads 

Hs = 4:Ec{n — UgY — Ej cos cj) + + 2/3eUrnis^(^i + (23) 


Here we have defined the quantized charge n on the superconducting island, which is proportional to the number of 
Cooper pairs, the offset charge Ug and the quantized flux (j) that threads the SQUID loop. The a{a^) are bosonic 
operators that act upon the resonator field, of frequency Ec stands for the charging energy of the superconducting 
island, while Ej = E^^^\ cos(7r^i/^o)| is the Josephson energy of the dc-SQUID loop. The coupling coefficient (3 is 
a function of the circuit capacitances, Ums is the root mean square voltage of the resonator, and e is the electron 
charge. 

One of the main problems that arise in using Cooper pair boxes as controlled quantum systems is the high sensitivity 
of their quantum state to noise on the offset charge Ug. By using large capacitive pads shunting the superconducting 
island, one can achieve insensitivity with respect to this noise. These Cooper pair boxes provided with large capacitive 
pads are known as transmon qubits [5]. In this parameter regime, the Hamiltonian of Eq. (23) models a transmon 
qubit. It can be rewritten in a diagonal basis for the part which is not coupled to the resonator. 


2 2 

Hs L^i \i) {i\ + uJrCi^a + gi^i-^i{\i) (i + I| + K + 1) (^|)(^ + <^^)5 

i=0 i=0 


( 24 ) 
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where we have included the first three levels i = {0,1, 2} of this device. Higher levels are usually excluded from the 
quantum dynamics, as being more energetic and hard to access. The first three levels can be approximately modeled 
as the first three levels of an anharmonic oscillator, with a relative anharmonicity factor of ar = (cj 2 — If ar 

is big enough, one can consider the transmon qubit device a two-level system [5]. 

The setups envisioned in this work use two variants of the transmon design. Each of the six devices sketched in 
Fig. la in the main text is made of three superconducting islands, connected via two SQUID loops, marked with 
different colors in the figure. The magnetic fluxes acting on these loops modify the quantum level structure of the 
system, and an appropriate manipulation of these fluxes in time allows for a fast modulation of the effective coupling 
to the resonator (see Eq. 24). This modulation does not affect the bare energy levels uji of the device [H [7] 
and gives the possibility of performing red and blue sidebands upon the coupled qubit-resonator system, which, in 
turns, allows to generate many-body interactions [8]. 

The setup considered in Fig. lb in the main text describes transmon qubits with improved design, the so called 
Xmon qubits [9]. In this case, the shape of the capacitance pads is designed in order to improve coherence times of 
the device, by reducing coupling to material defects and radiative losses. The interactions among these qubits are 
mediated by direct capacitive couplings, not making use of a common superconducting resonator. For this reason, 
the many-body interactions are reduced in this case to nearest-neighbor two-qubit gates. 


Many-body interactions and gauge invariance for different sectors 


In order to simulate the Hamiltonian associated with the non-Abelian model in Eq. (21), one can simulate stepwise 
the single many-body monomials that compose the total interaction. Each of them has the form of an A/'-qubit string 
term 


Hn = cr\cF^2 ' 


■ 


(25) 


where the set of superindices can assume values {i,j, ...k} G {x^y^z}. Up to a maximum of 2N local rotations, this 
general interaction can be mapped back to, e.g., erf erf • • • erf^. In what follows, we will show how to obtain this effective 
interaction in two different scenarios. In one case, one can use a set of N tunable-coupling transmon qubits coupled 
to a single microwave resonator. In the second scenario, a set of N charge-like qubits, e.g., transmon or Xmon qubits, 
are coupled to an additional ancillary qubit, for a total of + 1 qubits. 

We consider the sequence of collective and single qubit gates muni , as 


Us{t) 


^i(7r/4) Y.i<j g-d7r/4) 


and obtain the set of interactions 


Us{t) 

Us{t) =^ 9 tcrla^-crl ^ Au + 1, 

Us{t) = 4n, 

Us{t) AT = 4n - 2. 


(26) 


(27) 


The collective interaction 

by threading the two SQUID loops embedded in each device with fast-tunable magnetic fluxes, tuned at appropriate 


j in Eq. (26) can be realized within a set of N three-island transmon qubits. 


frequencies [8]. The general A^-body term of Eq. (25) can be simulated with this approach with a total of 2 collective 
gates, and a number of single qubit gates which is upper bounded by 2N -\-l. The Hamiltonian in Eq. (21) is composed 


of 16 monomials, 1 three-body term, 6 six-body terms, and 9 five-body terms. Therefore, the simulation for one digital 
step will amount to implement 32 collective gates and a number of single-qubit gates which is bounded by 184. 

One can implement the interaction in a set of capacitively-coupled charge-like qubits, by decomposing it into 
a set of C-phase and single-qubit gates. First, let us consider the sequence of gates in Eq. (26) between N qubits and 


an additional ancillary one, labeled by A. We consider also that the gates are rotated about the Y axis of an angle 


7r/4. The sequence of gates in Eq. (26) can then be rewritten in the following way 


UsaW = 




(28) 
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1.5 


FIG. 6. (Color online) Gauge invariance is preserved through digitization. The simulated gauge value converges 

to the expected invariant, 3/4 for a) and 11/4 for b), as the number of digital steps increases. The initial states |To) chosen 
in a) and b) are represented in the main text, inFig. 2 a) and c), respectively. The overlap |(T/|To)|^ witnesses the periodic 
oscillations of the initial state on the triangular plaquette for the first case, while it remains invariant in the second one. 


This sequence of gates is equivalent to the following interactions 

Us,a{^) =e^9tcTlat-a%^j^ ^ 

I7s,aW = + 

Us,A{t) AT = 4n, 

Us At) ^ 4^ _ 2 , 


(29) 


If the ancillary qubit is initialized in an eigenstate of cr^(cr^), |^a) = A)a)/^ a) ~ 

the effective dynamics on the A/'-qubit system would be equivalent to the interaction The interaction 

^j=i S can be realized as a sequence of N ZZ gates between the ancillary qubit and all the other ones, 
g2(7r/4)cr^ 1 cTj = ^ Each of these two-body gates can be realized as a sequence of two local Z- 

rotations and a C-phase gate between the qubits 1 and CPAA' Iii fact, defining 


g*0/2g-h0/2)cr^ _ 


one has that 



/ 1 0 0 0 \ 

/ 1 0 \ 

0 10 0 

( 0 = 

0 0 1 0 


( 0 0 0 / 

= e*' 

((/)/2)cr^crf ^ 


(30) 


(31) 


In this way, an arbitrary A^-body term can be encoded, up to a number of local rotations bounded by 2A^, with a 
sequence of 2N C-Phase gates and a number of single-qubit gates bounded by 6A^ + 1. Thus the quantum simulation of 


one digital step of the Hamiltonian in Eq. (21) will amount to realize 168 C-Phase gates and a number of single-qubit 
rotations which is bounded by 520. 

The C-Phase gate between two qubits can be performed by adiabatically changing the frequency of one qubit, 
bringing the level |11) close to the anharmonic excitated level |02), to induce state-dependent phase accumulations. 
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/ \ 

/ I I \ 



N / 

FIG. 7. (Color online) Scheme for an experiment involving two pure-gauge plaquettes, implemented in a lattice of capacitively 
coupled Xmon qubits. Multiple capacitive connections in superconducting pads with appropriate design can provide connectivity 
for larger lattices. 


Unwanted errors due to leakage to the |02) state are exponentially suppressed in the gate duration. For finite gate 
times errors can be minimized upon optimization in the Fourier spectrum of the control signal m- This optimization 
has been carried out for CZ gates in [12]. For a generic quantum simulation involving (7-Phase gates, one should 
define the set of angles to be used in the experiment, and optimize over the duration and Fourier spectrum of the 
gate, to avoid unwanted leakage. This has been done in m by increasing the gate time from 40 ns to 55 ns, not 
observing significant state leakages up to phases of 4 rads. Additional errors on the fidelity of the gate come from 
standard decoherence of the qubits, and control errors, see Ref. m for a discussion of the error budget. 


Finally, to complete the information given in the main text, we show in this section, in Fig. additional plots that 
complement Fig. 3 in the main text. The simulated gauge invariance is shown to converge for the singlet and triplet 
configurations on a single triangular plaquette. 


The proposed experiment can be extended to larger simulated lattices, which are implementable by considering 
setups of superconducting qubits with enhanced connectivity. The required connectivity can be achieved by using 
several microwave resonators, acting as non-local quantum buses, or by considering superconducting pads that can 
provide with multiple capacitive couplings, as schematized in Fig. For instance, we can envision, as shown in 
the main text, a ladder architecture that can support a singlet gauge invariant subspace, where non-Abelian string 
dynamics could be studied. 
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Scaling of errors for larger lattices 

It has been shown El that the scaling of the number of steps Ns needed to achieve a digital error e in a Suzuki- 
Lie-Trotter decomposition is almost linear in the simulated phase of the dynamics and sub-polinomial in e, according 
to the formula 


Ns < 


2m52*(m||iJ||t)i+V2fe 

gl/2/c 


(32) 


In the formula m is the number of terms that compose a generic Hamiltonian H = ^ fractal degree 

of approximation na, which we set to /c = 1. For one plaquette one has m = 16, 1 three-body term, 6 six-body terms, 
and 9 five-body terms, all present in the system with the energy J/2. Considering that the norm of an arbitrary 
product of P Pauli matrices is • • ’ CpW = 1, ^ {x,^, z}, one has, using a triangular inequality, that the 

norm of the SU{2) Hamiltonian for one plaquette is bounded by \ \H\ \ < \J\/S. As a consequence, the number of steps 
necessary to implement a dynamics with error e for the time t on a lattice of N plaquettes is bounded by 


Ns < 


2300N{N\J\t)^/^ 


(33) 


Notice that this formula gives a very loose bound which is not dependent on the underlying structure of the Hamil¬ 
tonian. Nevertheless, it shows that the scaling of the number of steps is efficiently solvable by a quantum simulator, 
not being exponential in the resources used. The exact scaling of the number of steps cannot be known (up to an 
exact solution of the problem, which is not addressable numerically for larger lattices), however it is useful to extract 
the dependence Ns ^ N^^^ from the upper-bound formula in Eq. (33). Assuming independent experimental errors 
on the single Trotter steps, one can relate the total number of steps Ns to the required experimental fidelity for the 
implementation of a single Hamiltonian term. We impose an accumulated experimental error of the same order 
of magnitude of the Trotter error ^ This enforces the error on the single step to be of the order ^ e/A^g. 
Besides these issues on error scaling and accumulation, the accesible simulation time has to be smaller than Ti and 
T 2 of the single qubits, defining a system of necessary conditions for the simulation of a lattice of N plaquettes 


< 


Ns 


^ 2300N{N\J\tf/^ 

- xA 


Ns 

U<Ti,T2. 


(34) 


Matter-field SU{2) Gauge Interactions 


In the main text we have focused on pure gauge dynamics with plaquette interactions, where usual perturbative 
approaches do not get energy scales which are big enough. In this section, we point out the basic elements needed for 
the simulation of the simplest SU (2) gauge model coupled to a matter field in (I + I)-dimensions. The basic ingredients 
and the validity of the different energy-scale regimes have been already discussed in the simulation of Abelian gauge 
models in trapped ions and superconducting circuits [num. This section can be seen as a generalization to non- 
Abelian setups. 

The ideal interaction Hamiltonian that we would like to simulate is given by 


Hint = E i + l)b{i + 1) + H.C. 

i 


= E^h*) 


CaL(*)c^i}(* + 1) 


E 

fiiy 


C'a/.cy (i)c^ij(i + l)cr^^ 


bi3{i + I) + H.C., 


(35) 


where ba{i) and Cap{i) are hard-core bosonic modes, with a G P ^ {L,R} and i G Z. 

In terms of these modes, the former Hamiltonian has two local (gauge) symmetries: 

• SU (2) gauge invariance around every site with generators: G^{i) = —I, i) + | Z^a /3 (0^S/3^/3(0+-C'“(*) * + l) 


• U{1) gauge invariance on the links with generator; G{i,i + 1) = J2a 


ciL(i)CaL(i) + + l)CaR{i + 1 ) 
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One way to obtain the previous Hamiltonian is to start with operators that respect one of the symmetries at the 
fundamental level and introduce an energy penalty for the combinations of operators that do not fulfill the second 
symmetry. In general, it is simpler to enforce the Abelian U{1) symmetry, since it only requires density-density 
interaction, and start with interactions that fulfill the non-Abelian symmetry from the beginning. 

The basic building blocks are given by a strong interaction that enforces the Abelian symmetry, i.e.. 


Ho=toJ2 


+ ^clp{i)c 

ap ('^) 


•n^[G(i,z + l)-7Vo]' 


(36) 


Here, the uj frequency of the bosonic modes splits the Hilbert space into subspaces with an equal number of excitations. 
Performing a rotating-wave approximation, interactions that conserve the number of particles are energetically allowed, 
while the term with energy scale Cl enforces Nq excitations per link. 

On top of the previous Hamiltonian of Eq. (36), a perturbative SU{2) invariant hopping term is applied 


^ = ^0 E I bii^CaLi^ + + y] + C^R{^<^f,aba{i) 


H.c. 


(37) 


Hence, in second-order perturbation theory, up to a renormalization of the bare frequencies uj of the modes, the 
effective Hamiltonian is given by 


= <1 \caL{i)4^{i + 1)1 + E + paA I bpii + 1) + H.c 

ia/3 pu ^ ) 


' y p 


(38) 


which reduces to the ideal Hamiltonian i^int, up to a density-density term that does not break any of the symmetries 
and it could be reduced or even eliminated, allowing a counter-term of the same type in the initial Hamiltonian Hq. 
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